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ABSTRACT 

We examine the relationship between star formation and AGN activity by constructing 
matched samples of local (0 < ^ < 0.6) radio-loud and radio-quiet AGN in the Herschel- 
ATLAS fields. Radio-loud AGN are classified as high-excitation and low-excitation radio 
galaxies (HERGs, LERGs) using their emission lines and WISE 22-(im luminosity. AGN ac¬ 
cretion and jet powers in these active galaxies are traced by [OIII] emission-line and radio lu¬ 
minosity, respectively. Star formation rates (SERs) and specific star formation rates (SSERs) 
were derived using Herschel 250-/im luminosity and stellar mass measurements from the 
SDSS—MPA-JHU catalogue. In the past, star formation studies of AGN have mostly focused 
on high-redshift sources to observe the thermal dust emission that peaks in the far-infrared, 
which limited the samples to powerful objects. However, with Herschel we can expand this to 
low redshifts. Our stacking analyses show that SERs and SSERs of both radio-loud and radio¬ 
quiet AGN increase with increasing AGN power but that radio-loud AGN tend to have lower 
SER. Additionally, radio-quiet AGN are found to have approximately an order of magnitude 
higher SSERs than radio-loud AGN for a given level of AGN power. The difference between 
the star formation properties of radio-loud and -quiet AGN is also seen in samples matched in 
stellar mass. 

Key words: galaxies: active —infrared:galaxies 


1 INTRODUCTION 


Since their discovery active galactic nuclei (AGN) have formed an 
important part of astrophysics research. Investigations of AGN are 
not only crucial in their own right, but also essential for galaxy 
formation and evolution studies. Accumulating observational data 
clearly show that the masses of black holes in massive galaxies are 
correlated with various properties of their hosts such as the galaxy 
luminosities (e.g. Kormendy & Richstone| 199^ [Marconi & Hunt 


2003[|Gultekin et al.|2009|l, the galaxy bu lge masses (e.g. |Magor- 

rian et al.|1998 McLure & Dunlop |2002|| and the velocity disper- 


(e.g'|Ferrarese & Merritt|2000[ Gebhardt et al.|2()00||Merritt| 


* E-mail: g. gurkan-uygun @ herts. ac. uk 


|& Ferrarese|2001| ). Furthermore, the anti-hierarchical evolution of 
AGN (e.g.|Miyaji et al.|2000||Hasinger et al.|2005[|Bongiomo et af] 
|2007[|Rigby et al.|2011| t, i.e. the fact that the space density of low- 
luminosity AGN peaks around z < I and that of high-luminosity 
AGN peaks around z ~ 2, is very similar to the cosmic downsiz¬ 
ing of star forming galaxies (e.g. |Cowie et al.|1996[ |Menci et~^ 
2008 [Fontanot et al.||2009|> an d spheroidal galaxies (e.g. |Cim^ 
et al. 2006[|Thomas et al.|2010| l. In addition to these, both the inte¬ 
grated cosmic star formation rate and the black hole accretion rate 
increase rapidly from z ~ 0 out to z ~ 2 ([Dickinson et al.|2003| 

Babic et af|2007l[Hopkins et al.|2007[[Willdns et al.|2008[[Rigby] 

et al.|201l| l. All these relationships indicate that the formation and 

growth of the black holes and their host galaxies are fundamentally 
intertwined. 
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Although these relationships are observed, there is not yet a 
clear understanding of how black holes grow, the link between the 
growth of black holes and their host galaxy properties, and what 
leads to these connections. To explain the observed co-evolution of 
black holes and their hosts, interactions (AGN feedback) between 
the black hole at the centre of a galaxy and the gas and dust that 


it contains have been invoked in theoretical models (e.g. Granato 
et al. 2004l|Springel et al. [2005] [Croton et al.|200^ [Hopkins et al. 


2006[|Sijacki et al.|2009[|Cen & Chisari|2011| l. Major mergers have 
been widely suggested as a triggering mechanism of AGN activity 
(e.g. [Kauffmann & Haehnelt|2000| [Hopkins et aH2005j >. Secular 
processes (disk instabilities, minor mergers, recycled gas from dy¬ 
ing stars, galaxy bars etc.) have also been discussed as a mechanism 
responsible for fuelling, in particular for low-luminosity AGN (e.g. 
[Johansson et al.|2()0^|Genzel et al.|2008l|Ciotti et al.|2010] l. 

Models often represent the effects of feedback in two ways, 
denoted ‘quasar mode’ and ‘radio mode’ (e.g. [Croton et al.|2006| l. 
In quasar mode the energy release occurs as winds with sub- 
relativistic outflows and wide opening angle driven by the radia¬ 
tive output of AGN. The generated radiation interacts with the gas 
and dust in the host galaxies and the resulting winds (either energy 
driven winds or momentum driven winds, see [Cattaneo et al.[ ( [2009[ ) 
for further information) can expel the gas from the galaxy. This can 
stop the accretion of matter onto a black hole and further quench 
the formation of stars (e.g. [Page et al.|2012[ |. However, studies of 


X-ray luminous AGN do not show any evidence for this (e.g. Har- 


[rison et al.|2012[ but s ee [Cano-Diaz et al.|20T^[Barger et al.|2014 
an d [Brusa et al.[2015[ l. 

In radio-mode feedback, the accretion of the matter does not 
lead to powerful radiative output; instead we see the production of 
highly energetic jets. The jets may play an important role in the 
fate of the host galaxy by heating up the cold gas and suppressing 
star formation (e.g. [Best et al.|2005[[Hardcastle et al.|2013[l. There 
is even an indication for this in radio-quiet quasars (e.g. Harrison| 
[et al.|2015[ l. In the most dramatic scenario, the jets can expel the 
molecular gas from the host galaxy and stop star formation. Radio 
mode feedback has been widely used in simulations as a mecha¬ 
nism to prevent the overproduction of stars (by shutting down the 
star formation) in massive galaxies in order to produce observed 
“red and dead” early-type galaxies (e.g. [Bower et al.|2006[[Werner| 
[et al.[[2013l [Chen et ^[2013b[ ). However, a positive radio mode 
feedback has also be en suggested ([Silk & Nusser|2010[ [Kalfount-"] 

[zou et al.[20I2[[2014l[Gaibler et al.[20I2| ). In this scenario, the radio 

jets drive shocks in the interstellar medium which enhance the star 
formation, something that has been observed for decades now (e.g. 


Bales & Rawlings|I990|[Best et al.| 1998 [[Jarvis et al. [200 1| [Inskip [ 
et al.|2005[ ). 

Different types of AGN might be able to provide feedback 
in multiple ways. It is known that low-excitation radio galaxies 
(LERGs) do not have radiatively efficient accretion, which removes 
the possibility of quasar-mode feedback. They are believed to be 
fuelled by advection-dominated accretion flows (e.g. [Narayan ^ 
[Yi|1995[ ADAFs) which can create an environment where the en¬ 
ergy release occurs kinetically by the radio jets. These jets are able 
to provide ‘radio-mode’ feedback. Radio-quiet active galaxies and 
radio-quiet quasars have typical AGN properties where we see the 
radiative output but no strong radio jets, and so we can only ex¬ 
pect to see ‘quasar mode’ feedback. On the other hand, radio-loud 
quasars and HERGs have radiative output produced by highly effi¬ 
cient accretion as well as kinetic energy release seen as strong radio 
jets. Therefore, both feedback mechanisms might be expected to be 
observed in these powerful objects. 


The anticipated relationship between AGN luminosity (or 
black hole accretion rate) and SFR has been investigated many 
times previously with mixed results (For recent reviews on black 
holes and galaxy evolution see [Alexander & Hick^[2012[ [Kor-[ 
[mendy & Ho[2013| and [Heckman & Best[2014[ l. For instance, some 
studies have found a slight correlation between these quantities 


(e.g. 

Lutz et al.|2008| Serjeant & Hatziminaoglou|2009|[Bonfleld[ 

[et al. 

2011 [[Silverman et al.|2009|[LaMassa et al.[2013|), some oth- 

ers found a strong correlation (e.g.[Hao et al.|2005|[Satyapal et al.[ 

12005||Netzer|2009| Shi et al.|20091|Chen et al.|2013a 

1 whereas oth- 

ers (e.g. Harrison et al.[2012 Rosario et al.[[2012 

Rovilos et al.[ 

[2012|[Mullaney et al.|2012| Bongiorno et al.|2012[l found a flat re- 


lationship (or no evidence for a correlation). On the other hand, the 
results of [Page et~^ ( [20I2| ) indicate a suppression of star formation 
due to AGN feedback (see [Caputi|2014[ for a recent review on the 
IR perspective of AGN). This wide variety of results is puzzling. 
It is important to note that although there are some overlaps be¬ 
tween samples, flelds and indicators used for SFR and AGN lumi¬ 
nosity, in general both the sample selection and the star-formation 
indicator used varies from study to study. These may cause biases 
in the conclusions derived. Thus, performing complete and coher¬ 
ent surveys of AGN, minimising the systematic uncertainties, and 
proposing revised models and testing them are crucial for future 
AGN research. Additionally, it has been pointed out that the differ¬ 
ent variability time scale of AGN activity and star formation can 
lead to these different results (e.g. [Hickox et al.|2014[ [Wild et ah] 
[2010[ [Hopkins & Hemquist|2009[ l. For this reason, instead of con¬ 
centrating on individual sources and instantaneous AGN activity, 
averaging the AGN luminosity over the populations should provide 
a clearer view of the relation between SFR and AGN activity (e.g. 
[Chen et al.[2013a[ l. 

Another interesting aspect is to search for differences in the 
relationship between the AGN activity and SFR for radio-loud and 
radio-quiet AGN. This has been investigated previously: [Dicken[ 
[et al.[ ( |20T2] ) found that at low redshifts radio-quiet AGN hosts have 
stars forming at higher rates than radio-loud counterparts while 
[Chen et ar] ( [2013b| ) observed two times more actively star-forming 
galaxies among radio-quiet AGN than radio-loud AGN for galaxies 
with stellar masses M>, > M©. These studies have provided 

important information on investigations of AGN—SFR relation for 
different types of AGN. However, it is still important to investi¬ 
gate similarities or differences between the relationship of AGN 
outputs (kinetic or radiative) and SFR using radio-quiet and radio- 
loud AGN samples matched in AGN power and stellar mass, and 
the possible reasons for these. 

In this paper we investigate the role of AGN activity in regulat¬ 
ing the host galaxy evolution as a function of the different types of 
AGN. In the pvQ-Herschel era, star-formation studies of radio-loud 
AGN concentrated on high-redshift sources in order to observe the 
thermal dust emission peaking in the far-infrared/sub-mm, which 
limited the samples to powerful objects (e.g. [Archibald et al.[2001[ 
[Reuland et al.|2004] ). However, it is now possible with Herschel to 
expand this to low redshifts ( e.g.[Hardcastle et al.[20I0[[Jarvis et al.] 
[2010[[Hardcastle et al.[20T3l[Virdee et al.[20I3[ l. Therefore, we are 

able to create matched samples of radio-loud and radio-quiet AGN 
with redshifts 0 < < 0.6 in terms of the relationship between 

their star formation properties and their AGN activity. 

The layout of this paper is as follows. A description of the 
sample and, the classiflcation of the AGN are given in Section 2. 
Our key results are given in Section 3, where the comparison of 
the stacking analysis between radio-loud and radio-quiet AGN as 
well as between HERGs and LERGs are presented. We also form a 


© 2011 RAS, MNRAS 000, [I]j^ 
































































































Radio-loud V5. radio-quiet: AGN activity and SFR 3 


sample of sources matched in their stellar mass in order to investi¬ 
gate the star formation properties of these sources by excluding the 
effect of mass. Additionally the relation between SFR and black 
hole accretion rate (Mbh) for radio-loud and radio-quiet AGN sam¬ 
ples is examined. In Section 4 we interpret our results. Section 5 
presents a summary of our results and conclusions. 

The cosmological parameters used throughout the paper are as 
follows: Qrn=0.3, Qa= 0.7 and Hq=10 km s“^ Mpc“^. 

2 DATA 

2.1 Sample and classification 

To construct our sample we selected galaxies from the seventh 
data release of the Sloan Digital Sky Survey (SDSS; |Abazajian| 
|et al.| 200 ^ catalogue with the value-added spectroscopic measure¬ 
ments produced by the group from the Max Planck Institute for As¬ 
trophysics, and the John Hopkins University (MPA- JHUf] It has 
31001 sources spanning a redshift range 0 < z < 0.7. The sam¬ 
ple does not include quasars and Type-I Seyferts because the AGN 
outshine the host galaxies for these objects which makes it difficult 
to study their host galaxy properties. The Herschel Astrophysical 
Terahertz Large Area Survey (H-ATLAS; |Eales et al.|2010| l North 
Galactic Pole (NGP) and the three equatorial Galaxy and Mass As¬ 
sembly (GAMA) fields ( [Driver et al.|2011| ) data were used to ob¬ 
tain far-IR fluxes of the sample galaxies. Stellar mass estimates are 
available for most of the sources in the catalogue (jKauffmann et al.j 
|2003b| ). 

[Best & Heckman| ( |2012| ) (BH12 hereafter) constructed a radio- 
loud AGN sample by combining the seventh data release of the 
SDSS sample with the National Radio Astronomy Observatory 


(NRAO) Very Large Array (VLA) Sky Survey (NVSS; ( jCondonj 
jet al.j 199^ ) and the Faint Images of the Radio Sky at Twenty cen- 


timetres (FIRST) survey ( 

Becker et al.|1995j) following the methods 

described by [Best et al.|( 

2005|) and Donoso et al.|(|2009|). We will 


briefiy summarize their method, further details are given by [Bestj 
jet al.j ( [2005] ) and jDonoso et al.|(2009j ). Firstly, each SDSS source 
was checked to see whether it has a NVSS counterpart: in the case 
of multiple-NVSS-component matches the integrated flux densi¬ 
ties were summed to obtain the fiux density of a radio source. If 
there was a single NVSS match, then FIRST counterparts of the 
source were checked. If a single FIRST component was matched, 
accepting or rejecting the match was decided based on FIRST prop¬ 
erties of these sources. If there were multiple FIRST components 
the source was accepted or rejected based on its NVSS properties. 
Since this classification is sensible and re-producible we use the 
classification of BH12 for our work here. 

We firstly cross-matched the initial sample of galaxies with 
the BH12 catalogue, and all objects that they classified as AGN 
form our ‘radio-loud’ subsample, which has 613 objects. The re¬ 
maining ‘radio-quiet’ AGN sample are classified using the modi¬ 
fied emission-line diagnostics given by jKewley et al.| ( [2006| ). This 
classification is shown in Figure Composite objects were sepa¬ 
rated from star-forming objects using a classification line given by 
jKauffmann et al.j ( [2003aj l. This classification line utilizes the [Nil] 
A6584 / [Ha] ratio so it is not shown in Fig.[^ 

The classification done by using optical emission line ratios 
and this process gave us 8035 star-forming objects, 1190 objects 


^ http://www.mpa-garching.mpg.de/SDSS/ 


Population type 

Counts 

Radio-loud AGN 

613 

Seyferts 

1190 

Composites 

2490 

Star-forming objects 

8035 

LINERs 

319 


Table 1. The number of sources in each population after the optical 
emission-line classification. 


classified as Seyferts, 2490 composite objects and 319 LINERs (Ta¬ 
ble [^. Emission lines were not detected for 17741 objects in the 
sample so they could not be classified in this way. This classifica¬ 
tion biases us towards massive, low redshift galaxies. The sample 
has a mean stellar mass, logio(M/M 0 ) = 9.55 and a mean red¬ 
shift = 0.14. After the classification using optical emission lines the 
mean stellar mass is logio(M/M 0 ) = 10.08 and the mean redshift 
equals 0.09 (Figure [^- bottom plot). The final sample has a red¬ 
shift range 0 < z < 0.56. We do not consider objects that are not 
classified using the BPT diagram (Fig. further in the paper, and 
we also discarded: star-forming objects, LINERs and the objects in 
the BH12 sample classified as star-forming sources from our AGN 
sample. The end result is that a combination of Seyferts and com¬ 
posite objects classified using optical emission lines formed our 
radio-quiet AGN sample, which has 3680 objects. 

The treatment of composite objects is discussed in detail in 
Section 3.1. For our analysis we use stellar mass estimations of the 
sources to be able to calculate their SSFRs. There are only 206 ob¬ 
jects in the sample do not have these measurements. These sources 
were excluded from the parts of the analysis that involve SSFRs. 
These sources are systematically at higher redshifts so the sample 
used in the parts of analysis that involve SSFRs will be biased. 

To be able to evaluate the properties of radio-loud AGN as 
a function of emission-line class, the radio galaxies in the sample 
were classified as high- and low-excitation radio galaxies (HERGs 
and LERGs). We initially used classification information provided 
by BH12. There are 191 sources which are not classified either as 
HERG or LERG by BH12. For these sources we used the WISE- 
based (details with regard to WISE data are given in Section 2.3) 
classification proposed by jGurkan et al.j \201A\ . Our final radio- 
loud AGN sample has 404 LERGs and 209 as HERGs. 

There are various diagnostics using radio and optical measure¬ 
ments to classify sources as radio-loud or radio-quiet AGN. We did 
not use any of these diagnostics to classify our objects as radio- 
loud or radio-quiet AGN. This is because our objects are Seyferts 
2s and radio galaxies and the estimation of optical luminosities 
for such objects, even when they are radiatively efficient, is diffi¬ 
cult. However, we checked whether our radio-loud and radio-quiet 
AGN classifications coincide with traditional classifications. In or¬ 
der to do this assessment, we used radio and emission-line measure¬ 
ments available to us. We calculated the ratio of radio luminosity 
(at 1.4 GHz) to an estimate of optical luminosity (at 5100 A). As 
our AGN are mostly Type II Seyferts, we do not know their intrin¬ 
sic optical luminosities and to derive these we used our bolometric 
luminosities {LhoU see Section 2.4) and the relation between Lboi 
and L 5100 given by jRunnoe et al.|p012j ). We have 1.4 GHz fluxes 
of the radio-loud sources from BH12. We cross-matched the cat¬ 
alogue of our radio-quiet AGN sample to the FIRST catalogue to 
obtain their 1.4 GHz fiuxes. This provided 131/1190 Seyfert detec¬ 
tions and 211/2490 detections of composites. For the sources that 
are not detected in FIRST we used the measured minimum fiux in 
the catalogue as an upper limit. We assumed a spectral index 0.8 
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Figure 1. In order to present the clear separation between Seyferts and Lin¬ 
ers we show here the [SII]/[Hq;] — [OIII]/[H^] diagnostic diagram for galax¬ 
ies in our sample, excluding sources from the radio-loud sample. A solid 
line represents the main AGN/SF division line summarized bylKewley et al.l 
( | 200 ^ . 


(Si, oc which is a typical expected value for radio sources 
(e.g. |Kellermann|19^|Hardcastle et al.|2010| ) and derived the 1.4 
GHz radio luminosities using FIRST fluxes of the sources. 

We exclude LERGs from the analysis because they do not 
have significant AGN-related optical continuum emission. The re¬ 
sults of this analysis are plotted in Fig. We see that there is a 
clear separation between the radio-loud and radio-quiet samples 
using this parameter, in the sense that the vast majority of ob¬ 
jects that we have classified as radio-loud following BH12 have 
r = logiQ(Lradio/Topticai) > 0.5, while almost none of the 
radio-quiet objects exceed this value. It is important to note that 
traditionally a threshold r > 1 is used to select radio-loud ob¬ 
jects (solid line on Fig. adopting this threshold would classify 
some of our radio sources as radio-quiet. However, (i) this is an 
essentially arbitrary dividing line, and (ii), given the systematic un¬ 
certainties in computing r for this sample introduced by the use 
of bolometric corrections, we cannot really be certain that the dis¬ 
crepancy between our observed threshold in r and the traditional 
one is significant. Reclassifying the very few radio-loud sources 
with r <C 0.5 as radio-quiet would not affect our analysis in any 
way. Given the broad consistency between the two radio-loudness 
estimates we choose to retain our original classification based on 
the BH12 analysis, which is easily reproducible by future workers. 
We also calculated two different radio loudness parameters which 
are commonly use d: the ratio of 5-GHz flux density to the flux den¬ 
sity at 2500 A (e.g. Stocke et al.|l992| ) an d the ratio of 5-GHz flux 
density to the flux density at 4400 A (e.g. Kellermann et al. 1989) . 
5-GHz radio flux densities were extrapolated assuming a spectral 
index 0.8. In the same way, we used the relations given by |Runno^ 
|et al.| ( [20T^ between Lboi and L 5100 , and Lboi and L 3000 to esti¬ 
mate fluxes at 5100 A and at 3000 A, respectively. Considering our 
sources are at low redshift these values will be approximately equal 
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Figure 2. Redshift (top panel) and stellar mass (lower panel) distribution 
of the whole sample before implementing the emission-line classification 
(black) and the distribution of the same properties of the sample after the 
classification process (grey). 


to the fluxes at 4400 A and 2500 A. By comparing these parameter 
estimations with the example shown in Fig. we found that there 
is an agreement between different radio-loudness ratios in terms 
of radio-quiet AGN having ratios lower than the given parameter 
value. For all the reasons mentioned above we did not use any of 
these parameters to classify our sources. All of our estimations are 
provided in the online table. We can also see from this figure that 
the redshift range is higher for radio-loud objects than radio-quiet 
sources. This is because there is not enough volume at low redshift 
coupled with the steep evolution of radio sources. Additionally, ra¬ 
dio sources’ host galaxies, which tend to be massive ellipticals, can 
be easily observed by SDSS to higher redshifts. 


2.2 Far-infrared data, star formation and specific star 
formation rates 

Herschel-AIhAS provides imaging data for ~ 600 square degrees 
using the Photo-detector Array Camera and Spectrometer (PACS; 
at 100-, 160-/im [rbar et al.||2010| [Poglitsch et al.||2010| | and the 
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Figure 3. The distribution of the radio-loudness parameter for the sample 
of radio-quiet AGN and HERGs used in this work. LERGs are not shown as 
they are not radiatively efficient AGN. We used an estimate of optical lumi¬ 
nosity at 5100 A and 1.4-GHz radio luminosity to calculate the parameter. 


Spectral and Photometric Imaging Receiver (SPIRE; at 250-, 350- 
and 500-/im, [Griffin et al.|20I0[|Pascale et al.|201I| Valiante et. al. 
in prep.). To derive the maximum-likelihood estimate of the flux 
densities at the positions of objects in the SPIRE bands the flux 
density from the point spread function (PSE)-convolved H-ATLAS 
images were measured for each source together with the error on 
the fluxes. Eurther details about flux measurements can be found in 
[Hardcastle et aL| { |20I0||2013] ). 

Ear-IR luminosity is widely used to measure star formation 
activity. As discussed by [Hardcastle et al.j ( |20I3j ) neither the in¬ 
tegrated far-IR luminosity nor the 250-/im luminosity is a com¬ 
pletely reliable estimator of SER because of the contribution by 
cold dust heated by old stellar populations ( [da Cunha et aL|20T0| 
[Smith et~n^[2012| >. However, we are secure using them as long 
as the far-IR emission is dominated by warm dust heated by star 
formation which, gives rise to dust temperatures ~25 K or more. 
We utilised 250-/im luminosity, as the star-formation indicator 
that is least affected by dust heated by the AGN (e.g. |Hatzimi-[ 
[naoglou et al.[[20I0| [Hardcastle et aL[[2010| [Barthel et al.[[2012| f 

Eor the sources with PACS and SPIRE detections we used a mod- 
ifled black-body spectrum for the far-IR spectral energy distribu¬ 
tion (SED; /j, oc —); we flxed the emissivity index ^ as 1.5 

ekT -1 

and allowed temperatures to vary to obtain the best fitting temper¬ 
atures, integrated luminosities (Lir) and rest-frame luminosities at 
250-/im (L 250 ) provided by the minimum values. To calculate 
the i^-corrections the same emissivity index (^=1.5) and the mean 
of the best-fitting temperatures for each class were used (See Ta¬ 
ble [^ for the temperatures). These corrections were included in the 
derivation of 250-/im luminosities. 

Table [^ shows the best-fitting temperatures for each type of 
AGN in the sample, together with uncertainties derived by the boot¬ 
strap method. The results show that the mean temperatures of the 
populations in the sample are around 25 K (Identical temperature 
estimates for star-forming galaxies were also reported by |Smith[ 
[et al.[[20T^ . This indicates that the far-IR measurements of our 
sources are not strongly affected by cold dust heated by the old 
stellar populations. It has been suggested that SERs derived from 
far-IR luminosity overestimate SERs of the sources with 
low M© yr“^ due to the contribution by old stars (e.g. 

[et al.[20I3] l. However, most of our sources in the sample have SERs 


SERbe- 

Clemens 


Population type 

Temperature (K) 

^r(K) 

Radio-loud AGN 

19.6 

1.06 

Radio-quiet AGN 

25.6 

0.37 

Composites 

27.2 

0.23 

Star-forming objects 

25.8 

0.12 


Table 2. The mean temperature values for different populations with errors 
on them derived by the bootstrap technique. 


above 1 M© yr“^ (Fig.[^ and any possible contribution by old stel¬ 
lar population should not be a problem for the SER estimates. 

Another source of contamination that we may observe at far- 
IR wavelengths is synchrotron emission from the jets of a radio- 
loud AGN. In order to check any possible synchrotron contamina¬ 
tion in the Herschel 250-/im band we cross matched our radio-loud 
AGN objects to the Giant Metre-wave Radio Telescope (GMRT) 
catalogue by following the method described by [Mauch et al.[ 
( [20T^ . This provided 325-MHz GMRT fluxes for 165 out of 613 
objects. The reason for this low success is that the GMRT survey 
only covers 36 per cent of the GAMA and NGP flelds, and has in¬ 
complete sky coverage and variable sensitivity. Assuming that this 
sub-sample is representative of the whole sample, spectral indices 
of the matched objects, derived using the NVSS and GMRT fluxes, 
were used to gain extrapolated fluxes at 250 /xm. Comparison of 
these with 250-/im far-IR fluxes showed that 100/165 objects have 
steep spectra and accordingly have much lower extrapolated flux 
than the measured far-IR flux at 250-/xm. 65 out of 165 sources 
have extrapolated fluxes higher than their far-IR fluxes at 250 /xm, 
but of these, almost all (59) are not detected at 250-/xm Herschel 
band. This demonstrates that the synchrotron contamination in the 
far-IR is not a serious issue for our sample. 

SER measurements for galaxies in the SDSS 7th data release 
catalogue, calculated using the H^ emission line, corrected for 
dust attenuation and flber aperture effects, can be found in the 
MPA-JHU data base ( [Brinchmann et al.[[2004| l. To determine the 
SFR/L 250 relation, we chose star-forming objects, classifled using 
emission lines (Section 2.1), and which have both L 250 and SFR 
estimates. We then used the median likelihood SFR estimates 
derived using the H^ emission line given in the catalogue to 
calibrate the relationship between SFR and L25c[3 The errors in 
SERs were estimated using the 16th and 84th percentiles. We 
then used a Markov-Chain Monte Carlo regression to obtain 
the relation between SFR and L 250 , considering the errors on 
both SERs and L 250 , and incorporating an intrinsic dispersion 
in the manner described by [Hardcastle et alT] ( [2010| ). The derived 
Bayesian estimate of the slope and intercept of the correlation are 

logio(-L 25 oAVHz-i) = 23.5(±0.01) + O.9O(±O.O2)logio(SFR/M0 
yr-i). 

The slope shown as the red line is very close to unity and 
the derived relation is in agreement with the results of [Hardcastle[ 
[et al.[ ( [20T3] ). The relationship between SFR and L 250 is presented 
in Figure [^ The relation between SFR and the integrated 8-1000 
/xm far-IR luminosity was also verifled and it is found to be in 


^ It would be desirable to use additional data to estimate the true SERs of 
galaxies taking into account the unobscured fraction that is not observed 
at far-IR bands. However, by calibrating our SER relation using SER de¬ 
termined using the Hq; emission line, we take into account the unobscured 
fraction in our estimates. 
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Figure 4. The distribution of the luminosity at 250-/im of the star-forming 
objects as a function of their SFRs. The best fit between the quantities is 
presented by a red line and icr errors of this ht are shown with grey dash 
lines. 


agreement with the results reported by |Kennicutt| l |1998| ) and 
[Murphy et al.| ( |20ir] ). 

For a further evaluation of L 250 as a star-formation indicator 
we compared the median-likelihood SFR estimates derived from 
Multi-wavelength Analysis of Galaxy Physical Properties (Mag- 
phys) SED fitting ( jda Cunha et al. 12008] [Smith et aL|2012'] l, where 
cold dust and unobscured far-IR emission are taken into account, 
and the SFR estimates from L 250 . This comparison showed that 
SFRs derived from these different methods are consistent with each 
other once the uncertainties are taken into account. In particular, we 
found that the slope of a regression line between the two SFR es¬ 
timates, and the mean ratio of SFRl 25 o to SVKMagphys, are both 
consistent with unity. 


2.3 Near- and mid-infrared data 

Mid-IR luminosity can be used as an AGN power indicator because 
emission due to accretion is obscured by dust and gas structure (if 
it is present) and re-radiated in the mid-IR. The WISE mission has 
observed the whole sky in four mid-IR bands (W1 [3.4-/im], W2 
[4.6-/im], FF3 [12-/im], WA [22-/im]) with an angular resolution 
of 6.1, 6.4, 6.5 and 12 arcsec, respectively. The WISE all-sky cata¬ 
logue was searched for all objects in our samples. This was done by 
searching the catalogue within 3 arcsec from the coordinates of the 
sample sources. 30416 sources of the sample before classification 
had WISE detections. The WISE measurements are given in Vega 


magnitudes so the magnitudes of the sources were converted into 
Jy using the standard WISE zero-point^ 

In order to compute the luminosity at 22-/im, 12-/im—22-/im, 
spectral indices (FI, a assuming a power law slope for the 
spectral energy distributions-SED) were calculated based on the 
WISE photometry for detected sources. K-corrections were not de¬ 
rived using SEDs because different components (old stellar popu¬ 
lation, torus and star formation) emit in the mid-IR. In order to be 
consistent we derived the index to be used for ^-corrections by 
taking the mean value of the indices of all sources. The mean index 
is 2.26 (cr=1.08). 

2.4 AGN power indicators 

The 22-/im and [OIII] A5007 fluxes and luminosities are available 
to us, and can both be used as AGN power indicators. In order 
to see from where the mid-IR emission originates we carried 
out Markov-Chain Monte Carlo regression analysis to obtain 
the relation between the L 250 (W Hz“^), as obtained from the 
temperature fitting, and the uLjy luminosity at 22-/im (L 22 erg 
s“^). Eor this analysis we used the star-forming objects that were 
detected in both Herschel 250-/im and WISE ll-jim bands. This 
provided the following relation: 

logio(F22/erg s“^) = -7.93 -f 1.18 (±0.04) log 10(±250/erg 

s-'). 

We then used this relation to derive the expected ll-fim lu¬ 
minosities due to star formation alone of all objects in the sample 
and the star-forming objects. In order to understand the effect of 
star formation on the ll-fim WISE band we subtracted these es¬ 
timations from the actual 22-/im luminosity measurements. These 
residuals were then plotted as a function of mid-IR luminosity. In 
Eigure we show the results of this analysis. In the top plot the 
relationship between logio(F 22 ) and logio(F 25 o) of the sources 
in our sample is seen, and we can see that there is a close rela¬ 
tionship between logio(F 22 ) and logio(F 25 o)- It is also apparent 
that radio-quiet and radio-loud sources have similar distributions to 
star-forming objects. The bottom plots present the residual 22-/im 
luminosities of star-forming objects (left) and sources of the radio- 
loud and radio-quiet AGN sample (right) as function of L 22 (erg 
s“^). The comparison of these clearly shows that the mid-IR emis¬ 
sion is dominated by star formation even for the AGN. As expected, 
star-forming objects cluster around the residual luminosity of 0. 

We implemented the same process using F[oiii] for a compar¬ 
ison, obtaining the relation for pure star-forming objects as follows: 

logio( L[oiii]/erg s“^) = -10.93 + 1.16 (±0.05) logio(L25o/erg 

s-i). 

In the same way this relation was used to derive F[oiii] for 
all objects in the AGN sample and for the star-forming objects. In 
order to understand the effect of star formation on the [OIII] lumi¬ 
nosities we subtracted these estimates from the actual F[oiii] mea¬ 
surements. In Eigurej^the results of this analysis are shown. In the 
top panel we see that there is a correlation between logio(F[oiii]) 

^ The relation used for the conversion can be found at 
http://wise2.ipac.caltech.edu/docs/release/allsky/expsup/ 
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Figure 5. The top panel shows the relation between L 22 and L 250 for star-forming objects and our sample sources that are detected in WlSE-22-iim band. In 
the bottom the distribution of the residual luminosities of the star-forming objects (left) and AGN in our sample detected at 22 iim (right) as a function of L 22 
together with a histogram of the counts are seen. The dashed line on the top panel shows the best fit for star forming objects. Star-forming objects are plotted 
as red crosses, radio-quiet AGN as blue squares, HERGs as open black circles and LERGs as open black triangles. 


and logio(L 25 o). However, radio-quiet and some radio-loud objects 
clearly have higher emission-line luminosities in comparison with 
star-forming sources. In the bottom plots, a boost in Iogio(L[oiii]) 
of the sources stands out when we compare with the distribution of 
star-forming objects seen in the left panel, although we see some 
trend of accumulation of sources around a residual of 0. To de¬ 
termine the cause of this the density maps for composite objects 
and radio-quiet AGN were produced. We also calculated the mean 
value of the residuals for both samples together with the error in 
these using the bootstrap method, which all are shown in Fig[^ It is 


clear in this figure that the composites peak around zero value of the 
residuals (the mean of residuals is 0.19 with cr=0.01) whereas the 
radio-quiet AGN are distributed above zero (the mean of residuals 
is 0.84 with cr=0.02). 

These results shows that 22-/im luminosity cannot be used as 
an AGN power estimator for this sample, (a similar result was re¬ 
cently reported by |Rosario et al. 12013) 1 but that [OIII] emission is 
not dominated by star formation. Therefore, we used [OIII] lumi¬ 
nosity in the calculation of the AGN powers. In order to compute 
radiative power we used the following relation: Pagn = 3500 x 
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Population type AGN Power 

Radio-loud AGN LERGs Jet power 

HERGs Jet power+Radiative power 
Radio-quiet AGN Radiative power 

Table 3. Table shows AGN power definitions for different type of objects 
used in this work. As previously explained, Jet power was derived using 
extrapolated 151 MHz flux densities of both HERGs and LERGs. Radiative 
power was derived using [OIII] emission line luminosities for radio-quiet 
AGN sample and HERGs. 


^[oiii] given by|Heckman et al.|(|2004|) which takes into account 
an average value for dust extinction. We can see in both Fig. and 
[^bottom left panels that there is a scatter to positive residual values 
between L 250 and L 22 , and L[oiii]. At 22 /im star formation dom¬ 
inates the band. Conversely, [OIII] is predominantly due to AGN 
with only a small contribution from SR 

To derive the jet powers, we utilised the relation Q(W) = 
SxIO^^x/L^g^, whereLisi is the luminosity at 15 1 MHz, in units 


, given by 


Willott et al. 


1999} and / is 10®/^ 


^151 

of lO'^^WHz-Sr _ 

(see |Hardcastle et al.|2007| ). We derived the fluxes at 151 MHz us¬ 
ing available fluxes at L4-GHz assuming the spectral inde^j^to be 
0.8. The jet power provides the best estimate of the AGN power for 
LERGs. It is worth noting that in principle the environment should 
be taken into account carefully in the conversion between radio lu¬ 
minosity and jet power (e.g. |Hardcastle & Krause|20T4] ). However, 
we cannot do this with the data available to us. For HERGs we 
combined the calculated jet powers and radiative powers to obtain 
the AGN power as discussed by [Mingo et al.| ( [20T^ . AGN powers 
used in this work are summarised in Table|3l 

[Boroson & Greenj ( |1992| ), jShen & H^ ( |2014| ) and [Sun &| 
jShenj ( |2015| ) have found that for a given AGN continuum lumi¬ 
nosity radio-loud AGN have higher [OIII] magnitudes or fluxes 
than radio-quiet AGN as a result of having different Eddington ac¬ 
cretion rates. Most of our radio-loud AGN contain LERGs which 
have different accretion mode than typical AGN. With regard to 
the SFR-AGN relation discussed in the following section HERGs 
only dominate the third and fourth AGN power bins. If we take 
into account the bias mentioned above we would expect HERGs to 
be included in AGN power bins lower than their current bins. This 
would not make much difference to our results since they would 
still have lower SFRs than their radio-quiet counterparts consider¬ 
ing the trend of their SFRs as a function of AGN power. 


3 RESULTS 

3.1 Stacking SFR: Radio-loud and radio-quiet comparison 

In Figure]^ we show the distribution of SFRs of the star-forming 
objects, and radio-loud and radio-quiet AGN sample against their 
stellar masses. From the flgure it can also be seen that radio-quiet 
and radio-loud AGN have lower SFRs compared to star-forming 
objects. Radio-loud galaxies not detected at 5cr in Herschel 250- 
fim band are indicated with arrows. 

To compare the star formation properties of radio-loud and 
radio-quiet AGN as a function of AGN power we used the Herschel 
250-/im band fluxes. However, the majority of the sources are not 
detected at the 5cr level in the Herschel 250-/im band. Therefore, 

^ In order to be consistent we use the same spectral index for all galaxies 
in the sample. However, we also used the GMRT estimates where available 
as a sanity check and found that this gave the same results. 


we stacked SFRs, SSFRs and stellar masses of the sources in 5 
AGN power bins. 

We determined the mean values of the SFRs and SSFRs of 
the samples in all bins. While calculating SFRs we first found the 
mean value of L 250 in each bin and then calculated the mean SFR 
in the corresponding bin by using the relationship derived in Sec¬ 
tion 2.2. This allows us to not be biased against sources that are 
weak or not formally detected. For the derivation of SSFRs, we 
divided the sum of L 250 by the sum of the stellar mass in individ¬ 
ual bins. This gave us the mean SSFR in the corresponding bins. 
The errors on these quantities were calculated using the bootstrap 
technique which enables us to derive the errors empirically, mak¬ 
ing no assumptions about the distribution of the luminosity. The 68 
percentile confidence intervals derived from the bootstrap samples 
were used as errors on SFRs. Errors on SSFRs were calculated by 
combining the uncertainties throughout the bootstrap analysis for 
SFRs and stellar masses in quadrature. All these estimates for both 
radio-loud and radio-quiet samples can be found in Table 

It can clearly be seen in Figure that SFRs and SSFRs in¬ 
crease with increasing AGN power for both radio-loud and radio¬ 
quiet samples. We also show the mean stellar masses for the bins 
to see the effect of mass on the observed relation between SFR and 
AGN power. Most of our sources in the radio-loud AGN sample are 
LERGs hosted by galaxies that are presumably massive ellipticals 
as indicated in (Fig. and [^. We see that radio-loud AGN re¬ 

side in more massive galaxies in comparison with their radio-quiet 
counterparts. Our results show that in general radio-quiet AGN of 
a given AGN power are found to have higher SFR and SSFRs than 
radio-loud AGN. The difference between the SSFRs of radio-loud 
and radio-quiet AGN goes up to an order of magnitude. 

Stacked measurements in confused images can be biased by 
the presence of correlated sources because the large PSF can in¬ 
clude flux from nearby sources. Several methods have been pro¬ 
posed to account for this bias, including the flux measurements in 
GAMA apertures by [Bourne et al.|p012| ). This method explicitly 
deblends confused sources and divides the blended flux between 
them using PSF information, so that average flux is conserved and 
not counted multiple times in the stack. We checked for the ef¬ 
fects of clustering in our stacks by comparing to the results ob¬ 
tained from average flux measurements in the catalogue given by 
[Bourne et al.[ ( [20I2] ), for the 1758 objects overlapping between the 
two samples. Average fluxes for the matched sample are slightly 
higher using the deblended apertures method compared with our 
stacking method, which indicates that our method is not biased by 
the effects of clustering and the results are robust. The reason for 
the slightly higher fluxes using the catalogue from Bourne et al. is 
likely because they account for extended flux outside of the cen¬ 
tral beam. These results are consistent with the trends in our full 
stacked samples. 


3.1.1 Ejfects of the intrinsic correlation between SFR and 
redshift, and the Malmquist bias 

It has been shown that SSFRs of galaxies evolve with redshift (e.g. 
[Oliver et al.[[2010| and references therein). The effect of this in¬ 
trinsic correlation found between SSFR and z should be taken into 
account. In order to check if our results are affected by this rela¬ 
tion we tried deriving the relation between SSFRs and z using our 
star-forming objects but found an unrealistically strong correlation 
between SFR and z. The reason for this is that the sample used 
here is biased in the sense of having star-forming objects that are 
selected using optical emission lines. This causes us to only se- 
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Figure 6. The top graph shows the relation between L[oiii] and L 250 for star-forming objects and our sample sources that are significantly (Sa) detected 
in their [OIII] emission-lines. In the bottom the distribution of the residual luminosities of the star-forming objects (left) and detected AGN in our sample 
(right) as a function of T^qiii] together with a histogram of the counts are seen. The dashed line on the top panel shows the best fit for star forming objects. 
Star-forming objects are plotted as red crosses, radio-quiet AGN as blue squares, HERGs as open black circles and LERGs as open black triangles. 


lect objects with high star formation at high redshifts. We are not 
able to match our samples in redshift and stellar mass because of 
the small size of the radio-loud AGN sample. Instead we used the 
relation given by |OIiver et al.| ( [20T0]) where they showed that late- 
type galaxies have SSFR oc (1+J)^®. We then took this relation 
into account and re-derived the relationships of SFRs and SSFRs 
as a function of AGN power (Fig. [^. Comparison of this with 
our initial results (Fig. clearly shows that our results are not af¬ 
fected by the intrinsic correlation between z and SSFR because the 
trends we initially observed between SFR/SSFR and AGN power 


do not change. Although we assume the same SSFR evolution with 
redshift for all types of galaxies here, radio-loud AGN have been 
found in more evolved galaxies in comparison to radio-quiet AGN 
and star-forming objects in the local Universe. Accordingly a less 
strong evolution is expected for these sources. Therefore, by as¬ 
suming a stronger correlation between SSFR and z for radio-loud 
AGN we may overestimate this effect on these particular galaxies. 

We note that Malmquist bias does not affect our derived SFRs 
in the stacking analysis since we include all objects, including non¬ 
detections, in the stacking. The parent sample is of course flux- 


© 2011 RAS, MNRAS 000,[I]|^ 











































10 


GUrkan et al. 



[OIII] Luminosity (erg s ') [OIII] Luminosity (erg s ') 

Figure 7. The comparison of density maps of the residuals (/o^fioCL^om] observed / L[oiii] measured)) both for the composites (left) and radio-quiet AGN 
(right) is shown. The mean values of these calculated for both samples are also shown together with the errors on them derived using the bootstrap technique. 


limited in [OIII] and/or radio luminosity, giving a correlation be¬ 
tween AGN power and redshift, but this is included by taking the 
SSFR—z relation into account described above. 

As indicated above due to limited sample size we are not able 
to carry out the stacking analysis by binning in both redshift and 
AGN power. However, we implement a separate analysis in order 
to check the effect of universal SFR/SSFR evolution with redshift. 
For this we estimated AGN power/SFR ratios for HERGs, LERGs 
and radio-quiet AGN. In Eigure m we show the distribution of 
the AGN power/SER ratio for all populations in our sample. Due 
to a large scatter, we binned the data in 3 redshift bins and cal¬ 
culated the median of the AGN power/SER ratio for each bin for 
all populations. Errors on these values were derived by the boot¬ 
strapping. In the first 2 redshift bins radio-quiet AGN have lower 
AGN power/SER ratio than HERGs and LERGs. This indicates 
that at all redshifts radio-quiet AGN have higher SER compared to 
HERGs and LERGs. Somewhat surprisingly we find that HERGs 
and LERGs have similar AGN power/SER ratios in most of the red¬ 
shift bins. Therefore, we conclude that the difference between SER 
of radio-loud and radio-quiet AGN is not due to the redshift effect. 


3.1.2 Effect of the composite objects 


Composites are sources that have emission-line characteristics in¬ 
termediate between pure AGN and star-forming objects, have in¬ 
termediate mass and SER and lie on a typical Baldwin-Phillips- 
Terlevich diagram (BPT; [Baldwin et al. 11981) betw een pure star¬ 
forming objects and AGN (e.g. [Salim et al.||2007| [Trouille et al.[ 
[20TT] ). We implemented the stacking process described above ex¬ 
cluding the composite objects from the sample to see if our results 
were affected by the inclusion of these objects. A striking differ¬ 
ence would mean that while selecting the composites we mostly 
pick sources with high levels of star formation. However, this com¬ 
parison showed that there is no difference between the SERs and 
SSERs of the samples with and without the composite sources 
compare Eigs.[^and[^. Eurthermore, [Elbaz et al.[ ( |20TT] ) demon¬ 
strated that the expected SSER of star-forming objects is around 


2.5 X10 yr ^. In both Eiguresj^and 1^ the highest mean SS¬ 


ERs are lower than this value. This also tells us that the compos¬ 
ites are not dominated by star forming galaxies. Therefore, we kept 
the composites in the radio-quiet AGN sample for all analyses pre¬ 
sented in this paper. Corresponding measurements are given in Ta¬ 
bled] 

3.1.3 Stacking SER: Radio-loud and radio-quiet comparison 
using mass-matched samples 

One of the possible effects that may cause the correlation between 
SER and AGN power we see in our results is the mass of the host 
galaxy. In order to judge its infiuence and to eliminate it, we im¬ 
plemented the same analysis for sources matched in their stellar 
mass. We carried out mass matching using the method described 
by [Virdee et al.[ ( |20T3] ). In summary, we randomly discarded 1 per 
cent of the galaxies of the comparison sample (the radio-quiet AGN 
sample). We then ran a K-S test to compare the distributions of stel¬ 
lar masses for the main sample (the radio-loud AGN sample) and 
comparison sample (the radio-quiet AGN sample). Repeating this 
process N times, where N is the number of comparison sources, 
enables us to select the best reduced catalogue. Once we obtain the 
best matching catalogue, 1 per cent of the remaining sources in the 
comparison sample is discarded and the process is repeated. We re¬ 
peat this process until the null hypothesis probability returned by 
the K-S test exceeds 10 percent. This mass-matching process was 
implemented for the sources in each AGN power bin separately. Al¬ 
though this method provides us the best matched catalogue, we lose 
the majority of the original comparison sample. The loss is around 
96 per cent so that our final sample of radio-quiet AGN contains 98 
objects. Eigure shows the distribution of stellar mass for each 
sample before and after the matching process. 

The results of this analysis are presented in Eigure Al¬ 
though there seems to be a slight increase in SER of radio-quiet 
AGN with AGN power the steepness of increase softened to a more 
gentle slope by leaving out the effect of galaxy mass. 

We also show emission-line radio galaxies separately in this 
figure. LERGs, in all bins, have the lowest rate of star formation 
as expected, since LERGs have been found in large ellipticals with 
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Figure 8. The distributions of SFRs of the sources in our sample detected in 250 fim Herschel band and star-forming galaxies plotted as a function of their 
stellar masses. Star-forming objects are plotted as red crosses, radio-quiet AGN as blue squares and radio-loud AGN as navy hlled circles. Radio-loud galaxies 
not detected at Scr the Herschel 250-/rm band are indicated with arrows. We do not show these for radio-quiet objects in order to present the distribution of the 
sources clearly as the number of radio-quiet AGN is high. 


little star formation (e.g. |Best et al.|2005] l. For a given AGN power 
HERGs’ host galaxies have comparatively higher star formation 
rates than LERGs but less than host galaxies of radio-quiet AGN. 


3.2 Correlation between black hole accretion rate and star 
formation rate 

If SER and AGN power are coupled then a correlation should be ex¬ 
pected between black hole accretion rates (Mbh) and star formation 
rates (Msfr) of these samples. In order to assess this, Mbh were 
calculated using Pagn=t]Mbhc^ where Pagn is the AGN power, 
whose derivation was discussed in Section 2.4 (see also Table[^, 77 
is an efficiency factor, Mbh is the black hole accretion rate and c 
is the speed of light. We assumed the value of efficiency to be 0.1 
for all sources. However, it should be noted that since the efficiency 
depends on the nature of accretion disk and accretion flows, a range 
of efficiencies can be found for different types of AGN. Especially 
in convection-dominated accretion flow (CDAE; e.g. |Narayan et al. 
2000 1 , advection dominated accretion (ADAE; e.g. |Narayan & Yi 
19951 and adiabatic inflow-outflow solution (ADIOS; e.g. Bland- 
ford & Begelman|1999| ) accretion flow models only a small fraction 
of the matter contributes to the mass accretion rate at the black hole 
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because of turbulence and strong mass loss. Eor this reason, much 
lower radiative efficiencies (77 0 . 1 ) are expected but the effi¬ 

ciency of the jet-generation process is not known. On the contrary, 
thin-disk accretion onto a black hole may lead to high efficiency 
factors (77 > 0.1, |Abramowicz et al.|1988HYu & Tremaine|200^ . 
The use of a single efficiency factor is therefore only an approxi¬ 
mation. 

A proxy of the rate of black hole growth (Mbh) as a function 
of star formation rate (a proxy of the rate of galaxy growth) is pre¬ 
sented in Eigure[^ where we see a strong relationship between the 
rate of matter accreted onto black hole and that of star formation 
for both types of active galaxy samples. 

We show the standard relation between bulge and black hole 
mass (e.g. IMarconi & Hunt|20^|Heckman et al.|2004| | with black 
solid line for a comparison to the relations we see between the rate 
of black hole growth and galaxy growth for our samples. When 
we compare our results with the observed standard correlation sug¬ 
gests that low-power AGN have black holes growing slower than 
expected to maintain the Mbh—MB uige relation (for example for a 
source with Msfr = 9x 10 “^ M© yr“^, the expected value of Mbh 
is around 10“^ M© yr“^) and black holes in high-power AGN are 
growing faster (for an object with Msfr = 2.5 x 10 ^ M© yr“^, the 
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Sample 

AGN power bins 

Mean 2 ; 

N 

Mean SFR 

AGN power bins 

Mean 2 ; 

N 

Mean SSFR 

Mean stellar mass 


erg s“^ 



M© yr-i 

erg s“^ 



(xlO-11 yr-i) 

(X 10+11 ]y[^) 

Radio-quiet AGN 

1.0e-r41-2.0e-r42 

0.03 

170 

(0.06) 

1.9e-r41-2.0e-r42 

0.03 

170 

4 QQ (0-50) 
(0.50) 

n 75 (0.03) 

(0.03) 


2.0e-r42-1.0e-r43 

0.06 

1099 

1 QQ (0-05) 
(0.05) 

2.0e-r42-2.0e-r43 

0.08 

1864 

Q 97 (0-28) 
(0.26) 

0 63 

(0.01) 


1.0e+43-1.9e+44 

0.11 

2082 

6 29 

(0.11) 

2.0e+43-2.0e+44 

0.13 

1333 

20 15 

n 87 


1.9e+44—9.8e+44 

0.14 

278 

1 1 '^7 

2.0e+44-1.0e+45 

0.15 

260 

25-69 

rj 09 (0.04) 

(0.04) 


9.8e+44-7.0e+45 

0.18 

50 

13 66 

(1.15) 

1.0e+45-7.0e+45 

0.18 

49 

30.86 

(0.08) 

Radio-loud AGN 

1.0e-r43-1.9e-r44 

0.04 

15 


2.0e-r43-2.0e-r44 

0.05 

16 

0.25 

2 64 

(0.29) 


1.9e-i-44-9.8e-i-44 

0.10 

57 

2 01 

2.0e-r44-1.0e-r45 

0.11 

60 

1 16 (0-36) 
(0.38) 

2 56 

(0.21) 


9.8e+44-7.0e+45 

0.25 

342 

019 (0.47) 

1.0e+45-7.0e+45 

0.22 

266 

1 70 

2 48 w-13) 

0.14 


7.0e+45-3.0e+46 

0.34 

178 

5 05 

(1-23) 

7.0e+45-6.0e+46 

0.26 

93 

9 59 (1-99) 
(1.05) 

9 75 (0.16) 

(0.18) 


3.0e-r46-6.5e-r47 

0.35 

20 

9c 90 (9-64) 

©Q 44) 

- 

- 

- 

- 

- 

HERGs 

9.8e-r44-7.0e+45 

0.30 

120 

4 85 

1.0e-r45-7.0e-r45 

0.23 

68 

1 43 (0-40) 

0.42 

3 05 ^0-41) 


7.0e-r45-3.0e-r46 

0.40 

75 

7 97 

7. Oe-t-45— 6 . Oe-i-46 

0.27 

22 

9 00 (1-^0) 
(1.65) 

9 91 (0.36) 

LERGs 

1.0e+43-1.9e+44 

0.05 

15 

0 66 

(0.51) 

2.0e+43-2.0e+44 

0.05 

16 

9 99 (0.25) 
(0.25) 

2 64 (O-Oi) 

(0.32) 


1.9e+44—9.8e+44 

0.10 

51 

1 QA (0-66) 
(0.70) 

2.0e+44-1.0e+45 

0.10 

53 

1 10 W-38) 

nio 

2 56 

0.23 


9.8e-r44-7.0e-r45 

0.23 

222 

9 99 (0.50) 

1.0e-r45-7.0e-r45 

0.21 

198 

1 67 

9 99 (0-11) 


7.0e-r45-3.0e-r46 

0.30 

103 

3 49 (1-59) 
(1.49) 

7. Oe-i-45— 6 . Oe-i-46 

0.25 

71 

1 91 (1-24) 
(1.29) 

2 70 (0-18) 

(0.19) 


3.0e+46-6.5e+47 

0.35 

13 

ll.OJ 

- 

- 

- 

- 

- 


Table 4. The results of the stacking analysis regarding our samples. In column 2 the chosen AGN power bins are shown. N indicates number of sources 
included in each bin. The mean measurements of SFRs in each bin with their errors calculated by the bootstrap technique are presented in column 5. New 
AGN power bins were defined for the mean measurements of SSFRs and stellar mass as the number of sources in the sample decreased by excluding sources 
that do not have stellar mass measurements. In column 6 the new AGN power bins are shown. The mean SSFRs and the mean stellar mass for each bin with 
their errors are shown in column 9 and 10, respectively. 


Sample 

AGN power bins 

Mean 2 ; 

N 

Mean SFR 

AGN power bins 

Mean 2 ; 

N 

Mean SSFR 

Mean stellar mass 


erg s“^ 



M© yr-i 

erg s“^ 



(xl0-“ yr-i) 

(xl 0 +“ M©) 

Radio-quiet AGN 

1.0e+41-2.0e+42 

0.02 

23 

-—( omr - 

1.9e+41-2.0e+42 

0.02 

23 

75710:23) 

0.24 

0 76 

0.10 


2.0e+42-1.0e+43 

0.05 

268 

1 00 

2.0e+42-2.0e+43 

0.06 

443 

9 95 (0.20) 
(0.21) 

0 76 

(0.03) 


1.0e+43-1.9e+44 

0.10 

657 

435 

(0.19) 

2.0e+43-2.0e+44 

0.11 

490 


0 73 W-03) 

(0.03) 


1.9e-i-44-9.8e-i-44 

0.14 

201 

7 48 (0-56) 
(0.56) 

2.0e-r44-1.0e-r45 

0.14 

191 


9 79 (0.05) 

(0.05) 


9.8e-r44-7.0e-r45 

0.17 

40 

11 80 

11.80 

1.0e-r45-7.0e-r45 

0.17 

40 

24 65 

0 72 

(0.08) 

Radio-loud AGN 

1.0e+43-1.9e+44 

0.05 

15 

0 66 

2.0e+43-2.0e+44 

0.05 

16 

9 99 (0.25) 

2 64 


1.9e+44—9.8e+44 

0.10 

57 

2 Qi (0.57) 
(0.61) 

2.0e+44-1.0e+45 

0.11 

60 

Y 26 (0-35) 

no5 

2 56 


9.8e-r44-7.0e+45 

0.25 

342 

919 (0.46) 

1.0e-r45-7.0e-r45 

0.22 

266 

1 70 W-27) 

2 48 


7.0e-r45-3.0e-r46 

0.34 

178 

5 05 

(1.29) 

7.0e+45—6.0e-i-46 

0.26 

93 

9 59 G-07) 
(1.12) 

9 7 c (O.I 6 ) 

(0.18) 


3.0e-r46-6.5e+47 

0.36 

20 

95 90 ( 11 . 10 ) 
ZJ..30 ( 10 . 59 ) 

- 

- 

- 

- 

- 


Table 5. The results of the stacking analysis of the sample where the composite objects are not included. In column 2 chosen AGN power bins are presented. 
N indicates number of sources included in each bin. The mean measurements of SFRs in each bin with their errors calculated by the bootstrap technique 
are shown in column 5. New AGN power bins were defined for the mean measurements of SSFRs and stellar mass as the number of sources in the sample 
decreased by excluding sources that do not have stellar mass measurements. In column 6 the new AGN power bins are shown. The mean SSFRs and the mean 
stellar mass for each bin with their errors are presented in column 9 and 10, respectively. 


expected Mbh is to be about 9x10 ^ M© yr ^). The implications 
of this result are discussed in Section 4.1. 

We also checked whether any possible contribution from star 
formation to the [OIII] luminosity can affect the relation that we 
obtained between SFR and Mbh- The red dashed line shows the 
expected relationship between accretion rate, derived by using the 
[OIII] luminosity, and Msfr for star-forming objects. The compar¬ 
ison of the red line with the results of stacking analysis for AGN 
demonstrates that the contribution from star formation to the [OIII] 
luminosity can only have an effect on sources in the lowest AGN 
power bin. 


Population type 

Mbh range 

SFR/Mbh 

Radio-quiet AGN 

lO-o - 2 x 10-0 

15000-1500 


2 x 10-0 - 2 x 10-1 

1500-100 

Radio-loud AGN 

lO-o - 2 x 10-1 

150-15 


2 x 10-1 - 2 xl 0 i 

15-2.0 


Table 7. SFR/Mbh ratios for radio-loud and radio-quiet AGN are given for 
different range of accretion rates. 
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Sample 

AGN power bins 
erg s“^ 

Mean 2 ; 

N 

Mean SFR 
M© yr-i 

Mean stellar mass 

(X 10+11 

Radio-quiet AGN 

2.0e-r42-1.5e-r44 

89 

0.12 

5 9 g (0-53) 
0.56 

9 no (0-11) 


1.5e-r44-4.5e-r45 

8 

0.18 

7 Qg (3.10) 

0 9 ^ (0.51) 

HERGs 

1.5e+44-4.5e+45 

59 

0.21 

0 00 (0.61) 

0 no (0.16) 


4.5e+45-5.3e+46 

36 

0.27 

6 45 G-70) 

(2.14) 

090 ( 0 - 19 ) 

LERGs 

2.0e+42-1.5e+44 

7 

0.04 

0 17 w-05) 

nSo 

1 12 

(0.12) 


1.5e-r44-4.5e-r45 

208 

0.17 

1 90 (0.22) 
(0.21) 

2 63 

0.08 


4.5e-r45-5.3e-r46 

105 

0.25 

3 60 

(0.97) 

o CO (0-11) 

(0.13) 


Table 6. The stacking analysis results of radio-loud (HERGs and LERGs) and radio-quiet AGN samples matched in their stellar masses. Similar to Table|^in 
column 2 chosen AGN power bins are shown. In column 3 the mean redshift measurements are given for each bin. N indicates number of sources included in 
each bin. The mean measurements of SFRs and stellar mass in each bin with their errors calculated by the bootstrap technique are presented in column 5 and 
6, respectively. We also show here emission-line radio galaxies included in each bin and their corresponding measurements as a reference. 


3.3 Stacking AGN power: HERG and LERG comparison 

The radio-loud AGN sample was investigated on its own to study 
the possible difference in the star formation properties of the 
sources as a function of emission-line type, discussed previously 
by [Hardcastle et al.| < |2010| |2013| ). Since we want to compare the 
SFRs and SSFRs for HERGs and LERGs, we defined new AGN 
power bins to include the similar numbers of sources in each bin 
for both classes. We present the results in Eigure In the top 
graph we can see that the change in the SERs is dependent on the 
AGN power. SERs increase, for both HERGs and LERGs, with in¬ 
creasing AGN power. The increase in SSERs for LERGs is more 
gentle in comparison to the increase of their SERs whereas we still 
see a clear increasing trend of SSERs for HERGs. In two out of 
three SER bins (and three out of four bins of SSER), LERGs have 
lower SERs/SSERs than HERGs, but the error bars are large. 

Since the size of the radio-loud AGN sample is small we only 
have four AGN bins for the stacking analysis of SER, and this is 
reduced to three bins for the estimates of SSERs when the galaxies 
with no stellar mass estimates were excluded. Therefore, a quantita¬ 
tive analysis was carried out to find the magnitude of any possible 
difference between the SERs of LERGs and HERGs. We divided 
the SERs of LERGs by the SERs of HERGs for each AGN power 
bin to calculate the SER ratio and the corresponding errors on these 
ratios. The mean ratio was then computed using these estimates. 
The mean ratio of SERlergs/SERhergs is 0.60zb0.32. This result 
indicates that within the errors there is at most a slight difference 
between the SERs of LERGs and HERGs of matched AGN power. 


4 DISCUSSION 

We have investigated the possible link between black-hole activ¬ 
ity and star formation by probing star formation properties of large 
matched samples as a function of AGN power. Results of this work 
indicate that AGN and star formation are coupled to some de¬ 
gree for radio-quiet and radio-loud AGN, and also for different 
emission-line types of radio-loud AGN. A contrast between the 
SERs and SSERs of different types of active galaxies stands out 
in our results. Host galaxies of radio-quiet AGN have more stars 
forming than their radio-loud counterparts. In the case of the dif¬ 
ferent emission-line classes, HERGs have higher SER/SSER than 
LERGs. 

The galaxy masses are different for these objects and the in- 
fiuence of this should be taken into account. Therefore, a stacking 
analysis has also been carried out for radio-loud and radio-quiet 


AGN samples that are matched in galaxy mass. Mass-matched sam¬ 
ples still indicate similar results; there is a slight increase in SER of 
radio-quiet AGN as the AGN power increases. This rate is higher in 
radio-quiet AGN than HERGs and LERGs, but an increase in SER 
and SSER is also seen for HERGs and LERGs. We found a wide 
range of Mbh which are correlated with the SERs of both samples. 
In the following section we interpret our results. 


4.1 Interpretation of the correlation between SFR and AGN 
power 

In the simplest interpretation the correlation found between the rate 
of star formation and black-hole activity could be taken to be ev¬ 
idence for the synchronised growth of black holes and their host 
galaxies (Eig. p^ . This would then suggest that the ex¬ 

isting cool gas supply on the host galaxy scale can feed the black 
hole in the centre of the galaxy at the same time as it allows for 
many more stars to form in the host galaxies of both radio-loud 
and radio-quiet active galaxies. In order to have this relationship, 
the available gas on host galaxy scales should be transported to the 
central regions by some physical mechanisms e.g. major mergers 
and secular processes. Major mergers are thought to be responsible 
for carrying a large amount of gas to the nuclear region of a galaxy 
(e.g. |Bournaud et al.|200^ [Hopkins et al.|2010] |, thereby trigger¬ 
ing powerful AGN and resulting in the co-evolution of black hole 
and its host galaxy in high accreting sources at high redshifts. Sec¬ 
ular processes (large galaxy bars, disc instabilities, minor mergers 
etc.) are considered to be sufficient to move some of the available 
gas to the inner regions of galaxies to feed black holes (e.g. ison| 
|erd]2^, especially in low accretion systems at low redshifts. 
Morphological studies of AGN host galaxies indicate that the cor¬ 
relations observed between the properties of black holes and that of 
their host galaxy are tight for ellipticals and bulge galaxies but not 
for disk galaxies with pseudo bulges (e.g. [Greene et al.|2008| ). The 
growth of ellipticals and bulge-dominated galaxies are commonly 
associated with mergers that can feed the most powerful quasars 
whereas low mass black holes have mostly been found in pseu¬ 
dobulges. The sources in our sample have low and moderate AGN 
powers at z < 0.6. We do not have complete and detailed mor¬ 
phological studies of all sources but our rough analysis (Section 
4.2) shows that most of the radio-loud AGN have elliptical hosts 
whereas radio-quiet AGN hosts tend to have more spiral character¬ 
istics. As mentioned above secular processes are often suggested to 
be responsible for fuelling disky (disk galaxies and pseudobulges), 
moderate luminosity galaxies in the low-redshift universe (e.g. [Ko-[ 
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Sample 

AGN power bins 

Mean 2 ; 

N 

Mean SFR 

AGN power bins 

Mean 2 ; 

N 

Mean SSFR 

Mean stellar mass 


erg s“^ 



M© yr-i 

erg s“^ 



(xlO -11 yr-i) 

(X 10+11 

LERGs 

7.0e+43-3e+44 

0.05 

22 

10 10 

COCO 

66 

00 

0 

7.0e+43-3.0e+44 

0.05 

22 

041 ^^ 

(0.19} 

269 ^^ 

(0.32} 


3.0e+44-1.5e+45 

0.13 

77 

1 9 q 0.43 

l.ZVo.42 

3.0e+44-l.le+45 

0.12 

54 

1 09 (0-37) 
n oi 

9 C9 (0.23} 


1.5e+45-6.5e+45 

0.24 

186 

2 71 

^0.52 

l.le+45-4.5e+45 

0.21 

149 

0 R7 

9 17 ( 0 . 13 } 


6.5e+45-2.3e+46 

0.30 

101 

9 C7I.46 

Z.O/^ 44 

4.5e+45-l.le+46 

0.25 

76 

1 RO 

9 09 (o.is} 


2.3e+46-7.0e+47 

0.34 

18 

13.56i:2i 

l.le+46-6.4e+47 

0.25 

37 

4 14 (2-54) 
(2.71} 

7 70 (0-23) 

HERGs 

3.0e+44-1.5e+45 

0.16 

23 


3.0e+44-l.le+45 

0.14 

10 

0 75 W-55) 

9 C9 ( 0 . 57 } 


1.5e+45-6.5e+45 

0.32 

96 


l.le+45-4.5e+45 

0.22 

51 

1 30 

3 05 


6.5e+45-2.3e+46 

0.39 

78 

8.62|;^| 

4.5e+45-l.le+46 

0.27 

26 

1 90 (0.71) 

3 24 


2.3e+46-7.0e+47 

0.40 

11 

35.25fi:2| 

l.le+46-6.4e+47 

0.28 

10 


2 41 (0-3i} 

(0.31) 


Table 8. Table shows the stacking analysis results for HERGs and LERGs. In column 2 the chosen AGN power bins are shown. N indicates number of sources 
included in each bin. The mean measurements of SFRs in each bin with their errors calculated by the bootstrap technique are presented in column 5. New 
AGN power bins were defined for the mean measurements of SSFRs and stellar mass as the number of sources in the sample decreased by excluding sources 
that do not have stellar mass measurements. In column 6 the new AGN power bins are shown. The mean SSFRs and the mean stellar mass for each bin with 
their errors are shown in column 9 and 10, respectively. 


[rmendy & Kennicutt|2004[ [Fisher & Drory|2011| l. Considering all 
these implications for the galaxies in our sample, it is likely that 
internal processes are in many cases responsible for the formation 
of stars and driving the gas into the nuclear regions to feed black 
holes in their centre. 


There is no clear answer to the question of when SFR corre¬ 
lates with AGN luminosity. However, the results of earlier studies 
indicate that a relation between AGN luminosity and SFR has been 
mainly observed for high-power sources (Lagn >10^^ erg s“^). 
Low-power (or quiescent galaxies) objects do not show such a re¬ 
lationship between these quantities in earlier works. For instance, 
|Diamond-Stanic & Rieke| \2012\ found no significant correlation 
between galaxy-wide SFR and Mbh in local Seyferts (the mean 
distance is around 22 Mpc). [Shao et al.| ( |2010| ), |Lutz et al.| ( |2010| ) 
and |Rosario et al.|p013| ) used hard X-ray selected local Swift-BAT 
galaxies (z < 0.3) and did not observe any correlation between 
Lagn and Lsf- Here in this work we also analyse a sample of low 
redshift AGN (z < 0.6). However, our work is different from that 
of |Rosario et al.| ( [2013| ) in that we take all components of differ¬ 
ent AGN types into account and calculate the total AGN power (by 
considering both mechanical and radiative output of AGN when it 
is present). We then derive average SFRs for our objects that are 
matched in total AGN power. Our analysis clearly shows that there 
is a relation between AGN power and SFR for mass matched sam¬ 
ples of high power and low power AGN (Fig.[^p^andp~5]). This 
is a manifestation of a connection between AGN power and SFR 
that holds for low- and moderate-power sources at low redshifts as 
well. As mentioned above, IShao et al.| ( [20T0l ), |Lutz et aL| ( [MTol ) 
and |Rosario et al.| ( [2013| ) have come to different conclusions: What 
leads to this difference? First of all, |Shao et al.| ( |2010] ) and |Rosarl^ 
|et al.| ( [20T3] ) used 60-/im far-IR luminosity, and |Lutz et al.| ( |2010] ) 
used 870-/im submilimeter luminosity as a SFR indicator, while 
we derive the SFRs using the 250-/im luminosity, which should be 
minimally contaminated by AGN, old-stellar populations and torus. 
Another important difference is that the host galaxies of X-ray se¬ 
lected AGN may be different from the host galaxies of mid-IR or 
optical selected AGN (e.g. |Hickox et al.|2009[ [Juneau et al.|2013| l. 
Furthermore, the hard X-ray luminosity comes directly from the 
accretion disk whereas the [OIII] luminosity is derived from the 
extended narrow-line region of AGN. For this reason, the hard X- 
ray is expected to show more variation on short time scales (years) 


compared to the [OIII] luminosity, introducing scatter into any re¬ 
lationship based on the power derived from X-rays. 

Another point to consider is the duty cycle of AGN. Some 
authors (e.g. [Hopkins et al.||2003] ) have attempted to explain the 
discrepancy between the observed Mbh-SFR relation and the ex¬ 
pected one in luminous objects using a simple duty cycle model 
in which it is assumed that the growth of black holes occurs at 
high accretion rates at a fixed Eddington ratio and time, and they 
are off (accreting at unobservably low rates) otherwise. However, 
Fig. shows clearly that black holes can grow slowly at lower 
accretion rates. Low-power AGN have black holes growing slower 
than expected from the Mbh—M suige relationship and black holes in 
high-power AGN are growing faster. There is actually a continuous 
smooth growth of black holes instead of an on and off phases and, 
therefore, our results cannot be explained with a simple duty cycle 
scenario. It is worth noting that with our available data we cannot 
examine whether the available gas reservoir for star formation is in 
the galaxy bulge or arms. Therefore, it would be interesting to im¬ 
plement the same analysis using SFRs of de-composed bulge and 
galaxy disk. 


Galaxy merger simulations have also predictions for the rela- 
tion between S FR and AGN activity. In these models (e.g. |Di Mat-| 


teo et al. 


2005) it has been suggested that the SFR/Mbh ratio ranges 


between 200—600 at the peak of star formation and this value is ex¬ 
pected to be around 1000 for low accretion rate objects (Mbh= 10“^ 
M© yr“^). Our results agree to some level with the predictions of 
this model. However, we have low and moderate luminosity galax¬ 
ies with spiral (mostly in the radio-quiet AGN sample) and ellipti- 


cal morphologies at low redshifts so many of the ellipticals in our 

sample may have ongoing mergers. [Ballantyne ( 

2008j) used star- 

burst disk models proposed by [Thompson et a . 

([2005 ) to study 


they found that SFR/Mbh should be around ~250 (for Mbh ~ 0.3 
M© yr“^) for a disk with Rout=l00 pc. This can only be observed 
for the whole galaxy of our radio-quiet AGN sample for black hole 
growth rates of Mbh = 2x10“^ — 2x10“^ M© yr“^. Recently, 
Thacker et al. |2014) examined the correlation between Mbh and 


SFR of all galaxies for various AGN feedback models using the 
time evolution of a merger simulation and they found that although 
most feedback models produce the observed M-a relation, there are 
distinct differences between the results of various models. Some of 
these models result in predictions that are in qualitative agreement 
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Figure 9. The results of stacking analysis of radio-loud and radio-quiet 
AGN. From bottom to top the graphs show: the comparison of SSFRs, stel¬ 
lar masses and SFRs of the samples stacked in their AGN powers. Since 
we included sources that do not have mass measurements in the stacking 
analysis of SFRs, this allowed us to have more sources and accordingly 
more bins for radio-loud AGN. The first two AGN power bins did not have 
enough radio-loud AGN but we still show the mean value of SFR corre¬ 
sponding to radio-quiet AGN. 


with our results reported here (see |Wurster & Thacker] ( |2013| ) for 
details of these models.) 


4.2 Probable reasons for the difference between SFR/SSFR 
of radio-loud and radio-quiet AGN 

As mentioned in the previous section, radio-quiet and radio-loud 
AGN present a moderately strong relationship between their SFRs 
and AGN powers. If the observed trends are examined in detail it 
can be seen that for radio-quiet AGN at the high power end there 
is a tendency of constant SFR with increasing AGN power. This 
may be due to AGN feedback that starts having a significant ef¬ 
fect on star formation when the gas supply reaches a certain level. 
In radio-loud AGN, black holes grow continuously. [Kauffmann &| 
[Heckman] ( |2009| ) also reached similar conclusions. They analysed 



24.40 


23.50 g 

i 

3 . 


22.60 


AGN power (erg/s) 


Figure 10. The results of stacking analysis of radio-loud and radio-quiet 
AGN where the intrinsic relation between 2 ; and SSFR given by jOliver et al.j 
12010) was taken into account. 


nearby galaxies selected from SDSS and found that the growth of 
black holes in local galaxies is present in two ways: if there is 
enough gas in the bulge the black hole can regulate itself; when 
the galaxy runs out of gas, the growth of black hole is regulated by 
the rate of mass loss due to evolved stars. 

Another question that arises from this work is: What could be 
the possible reason for the difference between the SFRs and SS¬ 
FRs of radio-loud AGN and their radio-quiet cousins? The obvi¬ 
ous difference between these AGN is the radio-emitting strong jets 
that we see in radio-loud AGN but not in their radio-quiet counter¬ 
parts. If we especially focus on the results we have from the mass- 
matched samples where we eliminate the effect of mass, the strik¬ 
ing difference in their SFRs may be suggesting that these strong 
jets are responsible for this difference, in the sense that there is the 
difference between the SFRs of radio-loud and radio-quiet AGN 
(lower/suppressed star formation in radio-loud AGN[^. 

Another reason for this difference may be due to the differ¬ 
ence between the morphology of galaxies in our mass sample. To 
check if we are comparing galaxies with similar morphologies, we 
searched the database of the Galaxy Zoo Project 2 (GZ2; |Willett| 
jet al.|201^ . This was done by scanning the GZ2 catalogue sources 


^ The relation between radiative power and SFR could also be examined 
as a further test. This can only be done using HERGs and radio-quiet AGN 
since only these galaxies have radiative outputs. However, the number of 
HERGs is not large enough in our sample to carry out a statistically signifi¬ 
cant analysis. It can be seen in Figurefldjthat we have only two AGN power 
bins with both radio-quiet AGN and HERGs. 
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Figure 11. The distribution of AGN power/SFR versus redshift for radio¬ 
quiet AGN, HERGs and LERGs. Black points show the median estimations 
of the AGN power/SFR ratio for all population for 3 redshift bins. Errors on 
these estimations were derived from the bootstrapping. 


with redshifj^z < 0.25. Although we did not have classifications 
for all of our sources, this process gave an idea of their morpholo¬ 
gies. We found that most radio-loud AGN are ellipticals but radio¬ 
quiet AGN have mixed morphologies (some of them are spirals 
and some are ellipticals). However, it should be pointed out that 
we cannot consider bulge and disk separately for radio-quiet AGN 
with our available data. 


4.3 The differences between SFR/SSFR in HERGs and 
LERGs 

We also had the advantage of exploring the star formation proper¬ 
ties of large samples of emission-line classified radio galaxies using 
SFRs estimated by 250-/im Herschel luminosity matched in their 
AGN powers (Fig. [^. It is apparent that there is a relationship 
between the star formation and black hole accretion rates. 

With regard to the comparison of SFRs and SSFRs for these 
galaxies, there is only very tentative evidence for a difference, in the 
contrast to the results of other studies (e.g. |Hardcastle et al.|2013[ 


Herbert et al.||2010[ |Baldi & Capetti|[2008| l. Previously [Hardcas^ 

tie et al. ( |2013| ) exWined star formation properties of HERGs and 
LERGs by binning in radio luminosity and found a clear difference 
between these galaxies. It is important to note that in the present 
work we have matched our samples in total AGN power rather than 
luminosity, to compare SERs of emission-line radio galaxies. Sepa¬ 
rating by radio luminosity (a proxy of the kinetic output) gives rise 
to a clear difference in the SERs of HERGs and LERGs whereas 
this is less apparent when we take into account both kinetic and ra¬ 
diative outputs of HERGs. We cannot rule out a model in which the 
AGN power—SER relation is the same for both radio-loud galaxy 
types, with the only difference being the way in which the accretion 
power manifests itself, though the HERG sample size in the current 
work is small. 

We also see a rise in SSERs of both samples with AGN power. 


® This redshift cut provides reliable classifications of galaxies, considering 
the SDSS image resolution. 



AGN power (erg/s) 


Figure 12. The result of stacking analysis of radio-loud and -quiet AGN 
where composite objects are not included. Similar to Figure|^from bottom 
to top the graphs show: the comparison of SSFRs, stellar masses and SFRs 
of the samples stacked in their AGN powers. 


When stellar mass is taken into account the rising trend of SSFRs 
of LERGs is not as sharp as we see for their SFRs. This may be 
expected, as we know that LERGs are hosted by massive red el¬ 
lipticals where most of the galaxy is dominated by an old stellar 
population, whereas HERGs tend to have lower mass and to be less 
evolved than LERGs. Other differences between these emission¬ 
line galaxies are the environment in which they are found (e.g. 
Hardcastl^|2004[ [Tasse et al. [2008 | l and the galaxy colours (e.g. 


Smolcic|2009[[Janssen et al.|20I2| |. If the HERG/LERG difference 


is understood as an accretion rate switch ( [Best & Heckman|2012[ 
[Mingo et al.||2014] | then we would expect LERGs to have lower 
accretion rates and/or more massive black holes, implying smaller 
amounts of cold gas and/or more massive (hence more evolved) 
host galaxies. 
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SDSS ID 

2 ; 

RA DEC 

F1uX 250 Err250 Spirst 

Flux ^01111 

Err^oiiii 

SFRh« 

P16sfrHc 

P84sfrHq 

Mass 

El^Mass 

Type 

588848899372744835 

0.07 

172.999 -0.598 

0.148 

0.006 

15.481 

3.027 

-0.310 

-1.050 

0.355 

11.177 

11.065 

11.264 

Comp 

588848899372744859 

0.11 

173.011 -0.534 

0.134 

0.006 

12.503 

2.158 

0.488 

0.119 

0.796 

11.009 

10.886 

11.114 

Comp 

588848899909746811 

0.12 

173.264 -0.154 

0.020 

0.006 

14.536 

2.664 

0.183 

-0.270 

0.626 

10.882 

10.773 

10.966 

Comp 


Table 9. Full tables of measurements of the sample are available online from the journal website. 




Figure 13. Stellar mass distribution of each AGN sample before (top plot) 
and after (bottom plot) the matching process. 

5 SUMMARY AND CONCLUSIONS 

In this work, we have explored star formation properties, derived 
using 250-/im Herschel luminosity, of radio-quiet and radio-loud 
AGN samples and investigated whether they depend on black hole 
activity. We have also compared these features of radio-quiet and 
radio-loud AGN samples matched in their stellar mass and red- 
shifts. The main results we have obtained are as follows. 

• Examination of radio-loud and radio-quiet AGN samples 
matched only in their AGN powers shows that for both AGN sam¬ 
ples the rate of star formation is increasing with increasing AGN 
activity. The same conclusion can be drawn for their relative galaxy 
growth rates (Fig.[^. 



AGN power (erg/s) 


Figure 14. The result of the SFR stacking analysis of radio-loud and radio¬ 
quiet AGN using mass-matched samples. The radio-loud objects are sepa¬ 
rated into emission-line class. 



Figure 15. The relation between the rate of black hole growth and the rate 
of galaxy growth is presented for our AGN samples. Assuming Msfr is 
proportional to bulge mass and Mbh is proportional to black hole mass, 
the black solid line represents the standard relation between the bulge and 
black hole mass. (e.g. [Marconi & Hunt|2()0^[Heckman et al.|2004) . Low- 
power AGN have black holes growing slower than expected value and black 
holes in high-power AGN are growing faster. The red dashed line represents 
the relationship between Msfr and Mbh that would be observed if star¬ 
forming objects were incorrectly classified as AGN. It can be seen that such 
contamination of the AGN sample cannot account for the observed trend. 
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AGN power (erg/s) 


Figure 16. The result of stacking analysis of radio-loud AGN sample. Dif¬ 
ferent AGN power bins were used in order to have the appropriate number 
of sources in each bin to be able to compare different emission classes. From 
bottom to top the graphs shows: the comparison of SSFRs, stellar masses 
and SFRs of the samples stacked in their AGN powers. 

• A comparison of the star formation properties regarding these 
samples indicates that the host galaxies of radio-quiet AGN are 
forming more stars for a given black hole activity than their radio- 
loud cousins. This difference in the level of star formation per unit 
stellar mass goes up to an order of magnitude (Fig.|^. 

• We also classified our radio-loud AGN sample in terms of 
emission-line type in order to compare their star formation prop¬ 
erties. Both LERGs and HERGs present the same trend of increas¬ 
ing SER with rising blaek-hole growth. In terms of their SSERs 
both types of galaxies present a similar trend; their SSERs are also 
increasing towards to higher AGN powers. They present almost a 
eonstant relative galaxy growth rate. There is only a marginal ten¬ 
dency for HERGs to have higher SERs and SSERs than LERGs 
when matched by AGN power (Eig.p^. 

• To aecount for the likely infiuenee of stellar mass we com¬ 
posed radio-loud and radio-quiet AGN samples matched in their 
galaxy masses. The same stacking analysis has been implemented 
for these samples. Our findings suggest that the black-hole growth 


and the star formation rate are coupled. The amount of stars form¬ 
ing for a given time increases with increasing AGN power for both 
AGN samples (Eig. [14}. 

• When we take into account the effect of intrinsic correlation 
between redshift and SER in our staeking analysis we have shown 
that the strong eorrelation between SER and AGN power for radio- 
loud and radio-quiet AGN sample does not vanish. Eurthermore, 
we have evaluated the effect of SER evolution with redshift. This 
analysis showed that the differenee between SERs of radio-loud 
and radio-quiet AGN is not due to this effect. 

• We have also assessed the relation between SER and Mbh and 
found that both radio-quiet and radio-loud AGN have a range of 
ratios of SER/Mbh- A comparison of our results with the observed 
correlation suggests that low-power AGN have black holes growing 
slower than expected from the Mbh—MB uige relationship and black 
holes in high-power AGN are growing faster. Our results cannot be 
explained with a simple model of duty cycle of AGN (Eig. ED. 

• Reasons behind the apparent difference in SER of radio-quiet 
AGN and HERGs and LERGs have been explored. This difference 
is still seen when comparing the samples of radio-loud and radio¬ 
quiet AGN matched in their stellar mass (Eig. E!}- There may be 
two possible reasons for this. Either the strong jets we observe in 
radio-loud AGN suppress star formation in these galaxies or the 
difference between the galaxy morphology of the AGN samples 
leads to this observed disparity, or some combination of the two. 
Nevertheless, a direet role for feedbaek from the radio jets is not 
ruled out. 
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